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VERSION WITH MARKINGS TO SHOW CHANGES MADE 

1. (Amended) A method for local and superficial [(probed volume and measuring 
area with dimensions on the order of one transport mean free path)] characterization of a turbid 
medium using the following parameters: 

1) the refractive index n of the turbid medium 

2) the absorption coefficient /x^of the turbid medium 

3) the reduced scattering coefficient ju/of the turbid medium 

4) the phase function parameter y -{l-g^W'^i) of the turbid media, where and are the first \ 
two moments of the Legendre polynomial development of the phase function /?, 

and comprising the steps of: 

- measuring the spatially-resolved reflectance R{p )oi the turbid medium (p being the source- 
detector distance) by any means [mean], comprising an illumination beam as a source and an 
optical detector, which, by using optional signal processing, which may involve filtering and de- 
convolution operation to correct for the non-zero area of either the illumination source or the 
detector, allows for the precise determination of the said spatially-resolved reflectance R{p), 

- mathematically processing R{p) to compute at least one of the said parameters: n, , /x/, y 
and/or the variations, in time and/or space, of at least one of the said parameters: An, A/x^, A /x/, 
Ay ,-whereby-an [the-said]""invef se problem", which consists in extracting the optical 
coefficients from the spatially resolved reflectance data is solved, and whereby a [the said] 
"direct problem" consists in computing the spatially resolved reflectance from the values of the 
optical coefficients n, /x^ , jit/, y involved in a model of propagation of the light in turbid medium 
and whereby the said "model" incorporates a Legendre polynomial development to the second 
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order of the said "phase function", and whereby the said "phase parameter" y is introduced in 
the computation as an independent parameter. 

2. The method of claim 1, wherein said spatially resolved reflectance is measured by 
a probe comprising at least one optical fiber carrying the light from the source unit to the turbid 
medium and at least one optical fiber collecting the reflected light and carrying it to detection 
unit, whereby the combination of a variety of emitting fibers and of receiving fibers yields a set 
of distances p at which the reflectance R{p) is measured. 

3. (Amended) The method according to claim [any one of claims] 1 [and 2], 
wherein said spatially resolved reflectance R{p) is measured for a set of values of p, by using a 
probe composed of optical fibers in any of the following configurations: 

one emitting optical fiber and a set of optical receiving fibers 

a set of optical emitting fibers and one optical receiving fiber 

a set of optical emitting fibers and a set of optical receiving fibers 
giving the spatially resolved reflectance R{p) at a variety of source-detector distances [distance] 
p and wherein the [. The] emitting and receiving fibers are [can be] arranged in one of the 
following [a variety of] configurations[. In some particular arrangement they can be positioned]: 

on a line, 

on crossed lines, 

on a circle, 

on an ellipse 

on crossed ellipses 
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on a disk, a rectangle, or any surface, as a dense and contiguous arrangement of fibers[, 
such as it can be obtained with a bundle of fibers or a multi-core fiber], 

on any pattern resulting from the combination of the above mentioned patterns. 

4. (Amended) The method according to claim [any one of claims 1, 2 and] 3, 
wherein said spatially resolved reflectance is measured by an optical and electronic micro-system 
comprising a collimated or focused beam as illuminating source and ID or 2D arrays of optical 
detectors[, arranged in a way similar to the one described in claim 3, and which, in some 
embodiment, can be a CCD or MOS camera]. 

5. (Amended) The method according to claim [any one of claims l-]4, wherein 
either [where] thiorobe [of claims 2 and 3] or the optical and electronic micro-system [of claim 
4] are put in ccmtact to the turbid medium. 

6. (Amended) The method of claim 1, wherein said spatially resolved reflectance 
is measured by a non-contact system, comprising at least one of the following combination of 
optical systems: 

a fixed optical system to irradiate the turbid medium with a collimated or focused beam 
forming the illuminating source and a fixed optical system comprising an imaging system 
fonning the image of the measured area of the turbid medium on a said "optical detector", which 
can be formed of ID or 2D array of optical detectors, whereby this second optical system can be 
identical to the first one and whereby the array of optical detectors can be either one of the 
following systems [system]: 
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■ a set of optical fibers [arranged in a variety of configurations sinailar to those described 
in claim 3], 

■ an optical and electronic micro-system (MOEM)[, which in some embodiment can be a 
CCD or MOS camera, with a functionality similar to the one described in claim 4], 

a fixed optical system for the collimated beam illuminating source and a scanned optical 
system for the said "optical detector", 

a scanning optical system for the collimated beam illuminating source and a fixed optical 
system for the said "optical detector", 

a scanning optical system for the collimated or focused beam used as an illuminating 
source and a scanning optical system for the said "optical detector". 

7. (Amended) The method according to claim [any one of claims] 1[. to 6.], 
wherein the absorption coefficient /x^ the reduced scattering coefficient fi/ and the phase function 
parameter y are determined by fitting the measured spatially-resolved reflectance R(p, /i/, ju^, y) 
to a set of discretized data obtained by using Monte Carlo simulations, or to interpolating 
functions giving a continuous approximation of the discretized data obtained by Monte Carlo 
simulations, and whereby said "Monte Carlo simulations" are based on a photon propagation 
model comprising a phase function approximated by a Legendre polynomial development 
limited to the second order[ In a preferred embodiment, cubic splines are used'as interpolating 
functions]. 

8. (Amended) The method according to claim [any one of claims 1 to] 7, wherein 
one or more of the following signal processing steps are [can be] performed [in addition to the 



- 14- 



BEVILACQUA et al 
Serial No. Unknown 



processing steps of any of the claims 1-6 whereby the processing of the spatially resolved 
reflectance data of claim 7 can be simplified and accelerated by]: 

fitting the measured reflectance R( p) by the function: 
rn^p"^ exp(m3p) 

to give [This fit gives] the values of the parameters m^ m^and m^, assuming that the expression 
R(p)= n\p"^- oxp(^m^p^ gives a smoothed expression of the spatially resolved reflectance R(p)[.]^ 

computing [Computing] the slopes — JR(p) and — (in R(p)) , or any mathematical 

dp ^ dp ^ ^ 

combinations of these two latter quantities and R(p), from analytical functions using the 

parameters m^ m^ m3, or by numerical procedures from the expression R(p)~ m^p"^ exp(^m^p^^ 

computing the values of at least one of the said parameters: n, /i^ , jit/, 7 and/or the 
variations, in time and/or space, of at least one of the said parameters: An, A/x^, A/x/, A7 from 
the relationship between the reflectance intensity R(p) and the slope of \nR(p) (denoted 
d p InR), determined at a fixed distance p comparable to the transport mean free path, whereby 
the computation is made from the data obtained by Monte Carlo simulations, provided that said 
"Monte Carlo simulations" are based on a photon propagation model comprising a phase 
function approximated by a Legendre polynomial development limited to the second order. 

9. (Amended) The method according to claim 7 [any one of claims 1 to 7], 
wherein the following signal processing steps are [can be] performed [in addition to the 
processing steps of any claim 1-6 whereby the processing of the spatially resolved reflectance 
data of claim 7 can be simplified and accelerated by]: 
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computing the reduced scattering coefficient /i/ and the phase function parameter 7 by 
using the following form of the reflectance: 

where the function A(p, 7, /i/j and B(ii^, fx/) depend also on the sources and detectors 
characteristics, and the refractive index of the sample, and comprising the steps of: 

computing [Computing] the slopes of the square root of the spatially resolved reflectance 



— V^Cp- Ai'/^vr y) = """(P^Mv'/)' which depends weakly on the absorption coefficient /x^for 
dp dp^ ^ 



d dA 

IP 0.3<p/x/<5, for at least two distances p, 

m 
W 

SI 



determining [Determining] the parameters /i/ and 7 by a polynomial interpolation of the 
data obtained by Monte Carlo simulations, whereby said "Monte Carlo simulations" are based on 
□ a photon propagation model comprising a phase function approximated by a Legendre 
polynomial development limited to the second order. 



10. (Amended) The method of claim [claims] 9, wherein the absorption coefficient 
/x^is determined by using the equation: 



where/ and h are continuous functions of the parameters p, 11/ and 7 [of] that can be 
approximated by a polynomial interpolation of the data obtained by Monte Carlo simulations, 
whereby said "Monte Carlo simulations" are based on a photon propagation model comprising a 
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phase function approximated by a Legendre polynomial development limited to the second 
order. 

11. (Amended) The method according to claim [any one of claims] 1 [to 7], 
wherein the difference Aix^=fx^ - /x^ between the absorption coefficient and a known value /x^is 

computed from the quantity yjRip, jX] , M«o ' 7 ) ~ V ^(P' l^l^l^a^Y) = B(a')' B(a^'), whereby the 

function B(a') of the albedo a' can be can be approximated by a polynomial interpolation of the 
data obtained by Monte Carlo simulations, whereby said "Monte Carlo simulations" are based on 
a photon propagation model comprising a phase function approximated by a Legendre 
polynomial development limited to the second order, and whereby the calculation can be done 
for a single y value, because the influence of the phase function parameter and 7 in B(a') are 
particularly weak. 

12. (Amended) The method according to claim [any one of claims] 1 [to 11], 
wherein the illuminating source is a broadband source or a white light source and the detector 
unit comprises a spectrograph or any wavelength analysis system to measure the wavelength 
dependence of at least one of the parameters (n, fi^fi/, y) • 

13. (Amended) The method according to claim [any one of claims] 4- [to4 2], 
wherein said turbid medium is a biological medium. 

14. (Amended) The method according to claim [any one of claims] 1 [to 13], 
wherein the measurement and processing is repeated at different locations [location] of the 
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sample [(multi-site measurements)], to build images of any one of the said parameters 

15. (Amended) An apparatus using tine nnethod of claim 1 [to 3, 5 or 7 to 14] 
for local and superficial characterization of a turbid medium, 

a) comprising a source, a detection unit, reference means, signal processing means, a probe 
comprising at least one optical fiber connecting said source unit to the turbid medium and at least 
one optical fiber connecting the turbid medium to the said detection unit, and reference means 

b) where the distance between the source and the detector is close to one transport mean free 
path^ 

16. (Amended) An apparatus using the method of claim 1[, 4, 5, or 7 to 14] 
for local and superficial characterization of a turbid medium, 

a) comprising an optical and electronic micro-system comprising at least one illuminating 
source, at least one detector, signal processing means and reference means, 

b) where the distance between the source and the detector is close to one transport 
mean free path 

17. (Amended) An apparatus-using the method of claim 1 [or 6 to 14] for 
local and superficial characterization of a turbid medium, 

a) comprising a collimated or focused beam used as an illuminating source, at least an optical 
detector for the detection unit, a fixed or scanning optical system for the illuminating source and 
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a fixed or scanning optical system for the said "optical detector", signal processing means and 
reference means, 

b) where the distance between the source and the detector is close to one transport mean free 
path 

18. (Amended) A system comprising three [Three] apparatus [similar to the 
apparatus] described in claim [any one of the claims] 15[, 16 and 17 ], and characterized by the 
fact that the distance between the collimated or focused optical beam used as illuminating source 
and the emitting point connected to an optical detector varies from 0.1 to 2mm. for application to 
biological media and to turbid media having a transport mean free path similar to biological 
media. 

19. (Amended) A test, where the control of the homogeneity of the sample over 
the probed area is performed with the apparatus of [any one of the claims] 15[, 16,17 and 18 ], 
which can be carried out according to the following [optional] procedure: in [In] the apparatus^ 
disposing two illuminating fibers [are disposed] symmetrically in regard to the collecting fibers; 
comparing [. If the sample is homogeneous,] the reflectance curves for each [curve should be 
identical with either] illuminating fiber to detect the f. Therefore,] heterogeneity of the 
investigated region or obstructions^ beneath the fibers : and, if [are detected by comparing the two 
curves. If] the two curves are sufficiently close, validating the measurement~[is.va]idated] and 
calculating the average of the two curves [is calculated]. 
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20. (Amended) A calibration and normalization procedure, which are carried out 
[optionally] with the apparatus of claim [any one of the claims] 15[, 16,17 and 18 ], whereby the 
following steps are performed: 

1) in [In] order to perform relative intensity measurements, the differences of transmitted 
intensity between each fiber for the apparatus of claim 15 [or the differences of signal intensity 
between each optical detector for any of the claims 16 and 17,] are determined by performing a 
measurement on a turbid phantom illuminated uniformly or a diffusing sphere of perfectly 
uniform properties;[. In this calibration procedure, the background light, measured with the light 
source turned off, must be subtracted from the signal. The obtained values are used to multiply 
the measurements given by each fiber and/or detector to correct the relative intensity 
measurements.] 

2) in [In] order to perform absolute intensity measurements, a calibration [can be] performed on 
a turbid medium of known optical properties, which can be realized according to any one of the 
following recipes [recipe] : 

a) solid or liquid turbid medium which properties have been measured by other techniques, 
or reported in the literature, 

b) water suspension of micro-spheres of known size distribution and refractive index[. 
Absorbing dye may be added to the suspension. In this case, the scattering properties are 
calculated from Mie theory, and the absorption coefficient is assumed to be equal to the water 
absorption coefficient, if no absorbing dye is added. If an absorbing dye is used, the absorption 
coefficient can be measured by a spectrophotometer, before mixing the solution with any 
scattering materials. 
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A Monte Carlo simulation is performed with the optical properties of the calibration sample. The 
simulation is then divided by the experimental reflectance performed on the calibration sample. 
The result, that must be independent of the source-detector separation, is defined as the 
calibration factor. Each new measurement can be multiplied by the calibration factor]. 
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